动态J积分及Rice J积分一致性比较与分析

首先回顾了动态J积分的发展进程,然后给出了两种常用动态J积分和Rice J积分的数学计算式,并对它们的一致性加以讨论,由此提出一种广义J积分的定义。然后通过算例,基于ANSYS软件,建立平面应变情形下纯I型中心裂纹模型,并利用编写的APDL程序计算了各个J积分,通过结果比较,进一步验证了在动态加载条件下,动态J积分比Rice J积分具有更好的路径守恒性,但仍可以看出它们的取值具有一定的一致性。最后,通过J积分间接法和裂尖位移外推法分别计算裂纹动态应力强度因子,所得结果吻合很好,验证了编写程序的正确性。     

边界元法分析薄形层合结构

薄形层合结构由于几何形状的特殊性,其力学分析是数值计算的难点.边界元法分析层合结构具有较大的优势,但对于薄形层合结构,边界源点和对边上的积分单元距离很近,边界积分方程中存在几乎奇异积分,常规的数值积分方法已经失效.文章引入一种半解析化方法,计算薄形层合结构边界元法中的几乎奇异积分,使边界元法能成功分析三维薄形层合结构的层间界面应力和各层内点力学参量.     

压电材料平面应力状态的直线裂纹问题一般解

研究了含直线裂纹系的压电材料平面应力问题单个裂纹和双裂纹问题的封闭解答表明，在裂纹尖端，应力、电场强度和电位移有１／２阶的奇异性并与前人结果比较了产生电场奇异性的物理因素     

汽轮机转子的温度场及热应力场计算

利用积分变换法导出了蒸汽温升率连续变化时汽轮机转子的温度场和热应力场计算公式。本公式可直接应用于转子的应力监测与控制。和差分法相比，解析法的数据处理速度快，能满足快速监测、及时跟踪与控制的要求。关键词：汽轮机转子；温度场；热应力场     

Stress Singularity in an Elastic Cylinder of Cylindrically Anisotropic Materials

The issue of stress singularity in an elastic cylinder of cylindrically anisotropic materials is examined in the context of generalized plane strain and generalized torsion. With a viewpoint that the singularity may be attributed to a conflicting definition of anisotropy at r=0, we study the problem through a compound cylinder in which the outer cylinder is cylindrically anisotropic and the core is transversely isotropic. By letting the radius of the core go to zero, the cylinder becomes one with the central axis showing no conflict in the radial and tangential directions. Closed-form solutions are derived for the cylinder under pressure, extension, torsion, rotation and a uniform temperature change. It is found that the stress is bounded everywhere, and singularity does not occur if the anisotropy at r=0 is defined appropriately. This revised version was published online in July 2006 with corrections to the Cover Date.     

Viscoelasticity: Fractal parameters studied on mammalian erythrocytes under shear stress

In this work, experimental determinations are carried out using a home-made device called an erythrodeformeter, which has been developed and constructed for rheological measurements on red blood cells subjected to definite fluid shear stress. A numerical method formulated on the basis of the fractal approximation for ordinary and fractionary Brownian motion¹ is proposed to evaluate the viscoelastic behavior of mammalian erythrocyte membranes. The diffraction pattern, which is circular when the mammalian erythrocyte membranes are at rest, becomes elliptical when the cells undergo shear stress. Photometric readings of light intensity variation along the major axis of the elliptical diffraction pattern are recorded during the creep and recovery process. These data series are used to calculate, fractal rheological parameters of self-affine Brownian motion on the erythrocytes, averaged over several millions of cells. Three different parameters over the time dependent process could be obtained, which are: correlation coefficient <C(t)>, correlation integral, andK ₂-entropy, and very different results were obtained.     

平面膜结构拓扑优化的有无复合体方法

将作者对桁架在应力约束下结构拓扑优化的有无复合体模型发展到平面膜结构在应力、位移约束下结构拓扑优化的建模与求解。同时提出了该模型的有效解法,获得了令人满意的数值结果。本文工作表明独立连续拓扑变量的提出对于结构拓扑优化的研究是有意义的。     

任意叠层球壳内的应力波传播

提出一种有效求解方法，求解了任意叠层球壳在动载荷作用下的应力波传播。首先，利用有限Ｈａｎｋｅｌ变换和Ｌａｐｌａｃｅ变换求解了每一单层的弹性动力学解。然后，利用叠层球壳的内、外边界条件和层与层之间的连接条件确定每一单层解中的待定常数。从而，我们可以得到应力波在叠层球壳内传播的精确解。最后，我们计算了在一个突加内压作用下的两层球壳（Ｓｔｅｅｌ／Ａｌ）和三层球壳（Ｓｔｅｅｌ／Ａｌｌｏｙ／Ａｌ）内的应力波传播。求解过程和计算结果表明文中提出的求解方法是简便可行的。     

弹性力学边界积分方程的一个发展

本文讨论二维弹性力学平面问题，独立于Ｒｉｚｚｏ型边界分方程，一类新型的边界积分方程，其边界场变量包含应力分量σｉｊｔｉｔｊ（其中ｔｉ是边界切向余弦）。该应力分量可直接用数值方法解边界积分方程求出，它比常规的边界元解提高一阶精度。文末的算例表明确定论的实用性和有效性。     

An analytical solution for the velocity and shear rate distribution of non-ideal Bingham fluids in concentric cylinder viscometers

An analytical solution is presented for the calculation of the flow field in a concentric cylinder viscometer of non-ideal Bingham-fluids, described by the Worrall-Tuliani rheological model. The obtained shear rate distribution is a function of the a priori unknown rheological parameters. It is shown that by applying an iterative procedure experimental data can be processed in order to obtain the proper shear rate correction and the four rheological parameters of the Worrall-Tuliani model as well as the yield surface radius. A comparison with Krieger's correction method is made. Rheometrical data for dense cohesive sediment suspensions have been reviewed in the light of this new method. For these suspensions velocity profiles over the gap are computed and the shear layer thicknesses were found to be comparable to visual observations. It can be concluded that at low rotation speeds the actually sheared layer is too narrow to fullfill the gap width requirement for granular suspensions and slip appears to be unavoidable, even when the material is sheared within itself. The only way to obtain meaningfull measurements in a concentric cylinder viscometer at low shear rates seems to be by increasing the radii of the viscometer. Some dimensioning criteria are presented.Notation A, B Integration constants - C Dimensionless rotation speed = µ/_y - c = 2µ - d = ₀ ²–2c_y - f() = (–₀)²+2c(–_y) - r Radius - r _b Bob radius - r _c Cup radius - r _y Yield radius - r ₀ Stationary surface radius - r Rotating Stationary radius - Y ₀ Shear rate parameter = /µ Greek letters Shear rate - = (r _y /r _b )²– 1 - µ Bingham viscosity - µ₀ Initial differential viscosity - µ µ₀-µ - Rotation speed - Angular velocity - Shear stress - _b Bob shear stress - _B Bingham stress - _y (True) yield stress - ₀ Stress parameter = _B +µ Y ₀ - _B - _y